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Abstract 

We establish in a cosmological model based on Lyra's geometry a relationship between a displacement 
vector field <p^, the Hubble parameter and the matter energy density p m (with a constant equation of state 
< uj m < 1) via an arbitrary function a(t). For a pressureless matter (ui m — 0) the effective equation of 
state oj e ff is completely determined by a(t). We subsequently investigate and find exact solutions in models 
yielding the ACDM and a power-law expansion. 
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I. INTRODUCTION 



Since the discovery of an accelerated expansion of the universe from astrophysical observations 
and the impossibility of explaining this acceleration with the Einstein gravitational field equa- 
tions of general relativity, R^ v — \g^ v R = T™, cosmologists think that, beside the standard matter 
there must be an unknown fluid with a negative pressure p DE < — |p DE , that drives the universe into 
acceleration. Such a component, dubbed "dark energy" was completely negligible in most of the 
past and will entirely dominate in the future. The origin and the exact nature of dark energy are 
still unknown. Subsequently, several interesting models have considered beside the standard mat- 
ter an alien dynamic dark energy component P-pT[ (for further reading see [1^, 12] and references 
therein). 

Attempting to solve geometrically the dark energy mechanism, a different approach considers the 
modification of the Einstein gravitational field equations with models involving only the standard 
matter. All the additional terms in the modified Einstein field equations beside the usual terms of 
general relativity constitute the so-called dark energy component. Such a modification can arise, 
either by extending the Einstein-Hilbert action to a more fundamental theory {fiR^R^^R^vap) 
theories of gravity constructed in pseudo-Riemannian geometry) or by modifying the Riemannian 
geometry. In the former case the simplest f(R) is that contains the cosmological constant A in the 
form f(R) = R — 2A. The cosmological constant thus introduces a constant contribution to the 
total energy- momentum tensor, which then has an effective equation of state uj e ff < — |, leading 
hence to an accelerated expansion. In this framework, the A dark energy is well-defined and reduces 
to a single number equivalent to a fluid with equation of state uj\ = p\/ pa = —1- Although the A 
dark energy is indeed consistent with all our current cosmological data and does not introduce new 
degrees of freedom, it encounters, however, the fine-tuning and cosmic coincidence problems by the 
fact that it is a non-dynamical quantum vacuum energy and its observed value is much less than 
the theoretical value computed from quantum field theory, p b s ~ 10~ 123 pl h ac . But several choices 
for f(R, R^iv, R^i/ap) have led to dynamic dark energy roughly compatible with observational data 



(for recent Review see [14, [15|] and references therein). Unfortunately, apart from f(R) = R — 2 A 
the modifications of the gravitational field equations from f(R, R^u, R^ap) theories include third 
or higher-order derivatives of metric components. Furthermore, higher-order terms in curvature 
invariants (such as R 2 , R^R^, the Kretschmann scalar R t j LVa ^R} 1,va ^ |16|], RDR or RO k R) or 
non-minimally coupled terms between scalar fields and geometry (such as <p 2 R) have to be added 
in an ad hoc way to the effective Lagrangian of gravitational field when quantum corrections are 
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considered. 

Another kind of modification of the Einstein gravitational field equations, mentioned above, is 



that constructed within Lyra's geometry [17], particularly in the so-called normal gauge. Lyra's 



geometry |17|-pi|1, like Weyl's geometry [22], is a modified Riemannian geometry proposed in order 
to unify gravitation and electromagnetism into a single spacetime geometry. However, in the 
framework of general relativity, only the connection in Lyra's and Riemannian geometries (contrary 
to Weyl's geometry) preserves the metric and the integrability of length transfers. Hence, Lyra's 
geometry becomes the most acceptable modified Riemannian geometry to describe the dynamic 
of the universe. Indeed, theories investigated in Lyra's geometry predict the same effects within 
observations limits, as far as the classical Solar System, as well as tests based on the linearized form 
of the field equations and are free of the Big-Bang singularity. They also solve the entropy and 



horizon problems, which beset the standard models based on Riemannian geometry [^, 23]. The 
modifications based on Lyra manifold are generally simpler than those built from f(R, Rfj, u , R^vap) 
theories, since they involve only first and second derivatives of metric components, and are basically 
due to the presence of an intrinsic geometrical vector field 4>^ (arising from the concept of a gauge) 
in the geometrically structureless manifold, without adding higher-order terms. On the other 
hand, the cosmological term {variable or constant), introduced in an ad hoc manner into the 
gravitational field equations has in Lyra's geometry an intrinsic geometrical origin arising from 
the displacement field [24-26]. In addition, the scalar-tensor theories of gravitation constructed 
within Lyra's geometry involve a scalar field and a tensor field that are all intrinsic to the geometry 

In this work we are to establish through an arbitrary function a(t) a connection between the 
Hubble parameter, the displacement field and the matter energy density, and express the effective 
equation of state in terms of a(t) and u) m . Considering a pressureless matter we also study a 
possible accelerated expansion of the universe, for a time-dependent and a constant functions a{t). 

The paper is organized as follows. In Section II, we describe a model, where a relationship 
between the Hubble parameter, the displacement field and the matter energy density is established, 
and we also investigate the ACDM model and Power-Law Expansion. Conclusions are given in 
Section III. 
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II. THE MODEL 



In Planckian units c = h = k 2 = 1, the Einstein gravitational field equations in normal-gauge 
Lyra manifold as obtained in [18] are: 



where G^ u = R^ v — \g^ v R is the Einstein tensor {R^ v and R being Ricci tensor and Ricci scalar, 
respectively), 

is an intrinsic geometrical tensor associated with the displacement vector field cj)^, whereas T^ v , 
defined by 

(3) 

is the energy-momentum tensor of a perfect fluid with energy p m and pressure p m = uj m p m (0 < 
Um < 1). The vector = (1,0,0,0), with u^u^ = —1, is the 4-velocity of the comoving observer. 
When the tensor Tuv vanishes, we recover the gravitational field equations of general relativity; so 
the modification of the Einstein's equations is solely due to the intrinsic geometrical vector field 

<?V 

In the flat Friedmann-Lemaitre- Roberts on- Walker background ds 2 = —dt 2 + a 2 (dr 2 + r 2 d0 2 + 
r 2 sin 2 9dip 2 ) , where a = a(t) is the scale factor of an expanding universe, the nonvanishing com- 
ponents of the Ricci tensor, and the Ricci scalar are given by 

R 00 = -3(H + H 2 ), Rij = (H + 3H 2 ) 9ij , R = 6(H + 2H 2 ), (4) 

where an overdot denotes differentiation with respect to the time coordinate t and H = a/a is the 
Hubble parameter. Throughout this paper we are to consider a time-like displacement vector 

<f>n = ((t)(t), 0,0,0). (5) 

Thus, the gravitational field equations (|l|) reduce to the Friedmann equations 

3H 2 + tj?{t) = p m , (6) 

- (2H + 3tf 2 ) +- A 4?{t) =u mPm , (7) 

while applying the Bianchi identities we find the total energy conservation equation 

p m + 3H(1 + oj m )p m = | (# + 3F</> 2 ) . (8) 



The matter energy interacts with the displacement field as the universe evolves. In [26] an inter- 
action was built between and the matter energy p m by encoding Eq. (JsJ) into two conservation 
equations. Here, we are to consider the first Friedmann equation (^) to find the connection between 
</>, H and p m . Assuming therefore the geometrical quantities H and <f) are both real, we can recast 
Eq. (|6|) in the form 

3(H+ l - ( j ) ){H- l -4>)=Pm (9) 

or compactly 

ZZ = \Z\ 2 = 1, (10) 
where the dimensionless quantities Z and Z are two c- numbers, functions of time t, defined by 

Z = J— (H + -</>) and Z = J— (H --(/>) with i = \/-L. (11) 

V Pm 2 V p m 2 

In terms of Z and its conjugate Z the Hubble parameter and the displacement field read 

H = (12) 
y 3 2 

and 

v V 3 2i V ' 

respectively. As for every c-number the general form of Z(t) can be parametrized as 

Z(t) = P(t)e ia{t \ (14) 

where a(t) and f3(t) are two arbitrary real functions of time. Hence, H and 4> become 

H(t) = ^(t)^Jt)co S [a(t)} (15) 

while 

4>{t) = -|/3(i)VMt)sm[«(t)]. (16) 

To reduce furthermore Eqs. (|l^) and (|l6|), we substitute them into (^) or use straightforwardly 
the identity ([H]) (ZZ = 1), and we find f3 2 (t) = 1. That leaves therefore a(t) as the only unknown 
function characterizing the c-number Z. Observationally, the Hubble parameter H and the energy 
density p m are positive; so taking (3{t) = +1 (from now on) the sign of H is determined by that 
of cos[a(t)], and we must have cos[a(i)] > whatever the sign of sin[a(t)] (eft can be positive or 



negative). Note also that, depending on sin[a(t)], the displacement field <fi could vanish yielding 
a decelerating expanding universe. But our goal is to explain the recent cosmic acceleration, thus 
we require sin[a(t)] ^ 0. 



Moreover, combination of Eqs. (15) and fllq ) relates H to <fi through the equation 



0(t) = 2tan[a(t)]H(t). 



(17) 



Thus, a(t) appears as playing the role of interaction between H, (ft and p, and will be called 
"interaction- function" from now on. 

On the other hand, in view of (|l7|), equation (^) transforms according to 



2H(t) + 3(1 - tan 2 : [a(t)])H 2 (t) = -0J m p m (t) 
while the energy conservation equation (||) becomes 

Pm(t) +3H[l + u m - (1 - w m ) tan 2 [a{t)]]p m (t) = 2a(t)t&n[a(t)]p m (t). 



(18) 



(19) 



Using the e- folding number x = In a, Eq. (19) may be written as a conservation equation of the 
matter energy p m with an 11 effective equation of state" Cj m {x) in the form 

dp m {x) 



dx 



+ 3 [l + u) m (x)] Pmix) = 0, 



where 



u> m {x) = oj m — (1 — tan 2 [a(x)] — - tan[a(x)] 



da{x) 



(20) 



(21) 



3 dx 

It is important to emphasize that, even if u> m = 0, the matter effective equation of state u) m (x) is 
nonzero, unless the interaction-function a(x) varies as 



a(x) = < 



mr, ui m = with n e Z, 



w. 



with x > - In (\k\), 



arcsin (fee 

arcsin (\J ui m + k e~ 3x ) , oj m ^ with < 



x > — In 

3 V u r, 

1, / k 
x > - m 

3 V 1 — UJr, 



if k < 0, 
if fc > 0. 



(22) 



k and k are constants of integration. 

The matter effective equation of state Q m (x) is not the effective equation of state of the system, 
Ueff, which by contrast is defined by 

2 H 



■1 



3 if 2 ' 



(23) 



By using ( jig ) and ( jig) knowing that /3(i) = +1, w e //(i) reduces to 

u eff(t) = u) m - (1 - U3 m ) tan 2 [a(t)]. (24) 

or equivalently, 

u e ff(%) = & m (x) + ~ tan[a(a)] d °[ X . (25) 
3 ax 

Which shows for tan[a(t)] ^ that the difference between u m and w e ff comes from the variation 
of the interaction-function. For a pressureless matter [oj m = 0), Cj m and uj e jf take the simple forms 



~ / \ 2r / m ^ da(x) 

Um{x) = — tan [a(x)J tan[a(x)j — - — (26) 

3 dx 

and 

u eff (t) = -tan 2 [a(t)], (27) 

respectively. So, when co m = both cu m and w e // may be completely determined by the interaction- 
function a(t) only. Assuming a positive H (then cos[a(t)] > 0) and any 0^0 acceleration occurs 
(for oj m = 0) when oj e ff(t) = — tan 2 [a(t)] < — |, i.e. for a(i) varying in the ranges 

-— < a(t) - 2n7r < — — or — < a(t) - 2nir < — with n e Z. (28) 
2 6 6 2 

In the following we are to consider a time-evolving and a constant interaction-functions and 

study then their consequences on the dynamics of the universe. 



A. ACDM model 

For oj m = and taking 

a(x) = arctan [(1 + Ae" 31 ')" 172 ] (29) 

with A a positive constant, leads as shown in |^] to the ACDM model and a constant displacement 
field 00) that plays the same role as the cosmological constant: A = |0g. It was found that |^6| 

Pm = Pcdm + 2A, 3H 2 = p CDM + A with p CDM = Pc DM a~ 3 , (30) 

where Pcdm is the present CDM energy density (that is, at t = to or x = 0). Notice that Eq. (ph 



can be computed straightforwardly by considering a constant displacement field 0(t) = 0o > in 



17|) and using (|i8[). One finds 

1 f 1 

- arccos ■ 

2 V! + Ce 



a(x) = - arccos ( - ; ) (31) 



with C being a constant of integration. After some algebras we will arrive at Eq. ( p9|) with the 
identity A = 2/(. 

Now we are to find explicitly the time-dependence of a(t) and quantities in the ACDM model. 



Using then ( \L7\ ) and (18) together with the identity A = 4^>q gives 



a(t) = arctan ( 1 = — ==- ) , (32) 

1 l + Ae^,'' 



where A is a constant of integration. Substitution of (32) into ( |i~7| ) and using the fact that H — > +oo 
at t = (yielding A = 1) enable one to obtain, 



ff(t) = x/A/s = A/A/Scoth (73A74t) (33) 

and 

a(t) = ai(l - e -^') 2/3 ev^* = 2 2 / 3 a, sinh 2 / 3 (y 7 ^*) (34) 
with cij a constant of integration. Also, cj e fy, tD m and p m are evaluated as 

" e//(t) = ~ U^ + J = " tanh (vW^t), (35) 



w m (t) = -1 + ^=- t = -1 + sech(V3At) (36) 

1 — I - 6 



and 



, 1 

p m (t) = 2 A — — 2 = Acosh (v3Ai)csch 2 (V3A/4i). (37) 

( e V3Ai _ i^j 

Eq. (|37| ) shows that p m {t) 2A when i becomes very large, and only the cosmological constant 
contributes to the matter energy p m (t), as seen from (|30|), 

On the other hand, Eqs. ( j30| ) and (|34|) allow one to compute the CDM energy p CDM and the 
constant of integration af. 



p CDM (t) = Acsch 2 ( v / 3A/4t) (38) 

and 

3 _ I CDM _ J CDM /on\ 

l ~ 4A -4(1-O0 dm) < ( 39 ) 
where Qq DM = p® dm /3Hq is the present fractional density of CDM. 



B. Power-law expansion 



Here, we are going to study the case ll co e ff(t) = u m (t)" by considering a nonzero constant 

interaction-function, a(t) = «o- This leads for oj m = to a power-law expansion (derived from 
(ID): 

H(t) = - f Kt~^ ~ OI a(t) OC ^/[3(l-ta^ 2 (ao))]_ (4Q) 

w 3[1 - tan 2 (a )] t 

Other cosmological quantities will be given by 

Ueff (t) = u m (t) = - tan 2 (a ), 0(t) = \ (41) 



. . 4 1 + tan 2 (a ) 1 

and ^ = 3 [l-tan 2 (aoF (42) 

Assuming -ff(i) and <j>{£) are all positive requires the conditions tan(«o) > and tan 2 (ao) < 1, 
what is translated by 

3"7T 7T 

— 7r < ao — 2n7r < — — or < ao — 2n7r < — with n e Z. (43) 



In addition, the cosmic age problem is alleviated when together with the constraints fl43|) ao a ls° 
satisfies 2/[3(l — tan 2 (ao))] > 1, that is, 

< ao — 2n7r < — — or — < ao — 2n7r < — with n e Z. (44) 

6 4 6 4 

Hence, as the ACDM model, the a^-model may satisfy the stellar bound on the age estimation 
[]34] |36|] and is free of the problem of the age of the universe [37]. The conditions ( f44| ) are exactly 
those yielding an accelerated expansion of the universe. 



III. CONCLUSIONS 



In this work we constructed within Lyra's geometry a cosmological model involving a matter en- 
ergy, with equation of state < oj m < 1, interacting with the displacement field. After establishing 
a connection between the Hubble parameter and the displacement field through an interaction- 
function, we showed that the conservation equation of the total energy in the universe (matter + 
displacement field) could be written as that of the matter energy with an effective equation of state. 
Having considered, furthermore, a pressureless matter we studied an accelerated expansion for a 
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time-dependent and a constant interaction-functions, leading to the ACDM model and a power-law 
expansion, respectively. 
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